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Abstract: Wingless-type (Wnt) signaling proteins participate in various cell developmental 
processes. A suppressive role of Wnt5a on keratinocyte growth has already been observed. 
However, the role of other Wnt proteins in proliferation and differentiation of keratinocytes 
remains unknown. Here, we investigated the effects of the Wnt ligand, Wnt3a, on proliferation 
and differentiation of keratinocytes. Keratinocytes from normal human skin were cultured 
and treated with recombinant Wnt3a alone or in combination with the inflammatory 
cytokine, tumor necrosis factor a (TNFa). Furthermore, using bioinformatics, we analyzed the 
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biochemical parameters, molecular evolution, and protein-protein interaction network for 
the Wnt family. Application of recombinant Wnt3a showed an anti-proliferative effect on 
keratinocytes in a dose-dependent manner. After treatment with TNFa, WntSa still 
demonstrated an anti-proliferative effect on human keratinocytes. Exogenous treatment of 
Wnt3a was unable to alter mRNA expression of differentiation markers of keratinocytes, 
whereas an altered expression was observed in TNFa-stimulated keratinocytes. In silico 
phylogenetic, biochemical, and protein-protein interaction analysis showed several close 
relationships among the family members of the Wnt family. Moreover, a close phylogenetic 
and biochemical similarity was observed between WntSa and WntSa. Finally, we proposed a 
hypothetical mechanism to illustrate how the WntSa protein may inhibit the process of 
proliferation in keratinocytes, which would be useful for future researchers. 

Keywords: keratinocyte; tumor necrosis factor; Wnt signaling 



1. Introduction 

The Wingless-Type MMTV (mouse mammary tumor virus) integration site family (Wnt) of 
signaling proteins participates in various developmental process, such as controlling cell proliferation, 
cell-fate determination, and differentiation during adult homeostasis [1]. Classically, on the basis of 
their capacity to inhibit phosphorylation of p-catenin by glycogen synthase kinase (GSK)-3p and its 
subsequent degradation, they are termed as canonical (e.g., Wntl, 3a, 8), whereas if they do not have 
any effect on P-catenin levels then they are termed as non-canonical (e.g., Wnt4, 5a, 11) [2,3]. 
The canonical Wnt ligand, exemplified by Wnt3a, binds to Frizzled (Fz) receptors and co-receptors, 
low density lipoprotein receptor-related protein-5 (LRP-5), or LRP-6, leading to the activation and 
nuclear translocation of P-catenin. In the nucleus, P-catenin acts as a transcription factor and governs 
key cellular processes, such as cell fate determination and stem cell maintenance [4]. By contrast, 
activation of non-canonical Wnt signaling by ligand like WntSa is involved with frizzled (Fz) receptors 
in conjunction with alternate co-receptors, including RORl/2 or Ryk, causing P-catenin-independent 
changes, such as protein kinase C activation and cytoskeletal rearrangements [5]. However, participation 
of Wnt proteins in the activation of signaling pathways is not as simple as it is presumed. Recent 
reports on the involvement of WntSa [6] and Wntl 1 [7] in the activation of P-catenin signaling suggest 
activation of multiple pathways, depending on which Fz receptors are expressed on the cell surface [3]. 
In addition, various antagonists of Wnt signaling, such as secreted Frizzled-related proteins (sFRP) [8], 
Dickkopf (Dkk) proteins [9], Wnt inhibitory factor (WIF)-l [10], create an extremely complex system of 
action for Wnt signaling. 

Psoriasis is a chronic debilitating disease of the skin characterized by sharply demarcated scaly red 
plaques commonly located on the extensor surfaces of the skin. The most characteristic feature of 
psoriasis is chronic inflammation, abnormal proliferation, and differentiation of keratinocytes [11]. 
The psoriatic epidermis exhibits aberrant expression of antigens associated with hyperproliferation, 
such as the heterodimer keratin 6-keratin 16 and heat-shock proteins [12]. Moreover, there is some 
evidence that points to immune activation with lymphocyte infiltration in the psoriatic epidermis and 
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secretion of inflammatory cytokines, such as interleukin (IL)-17, TNFa, and interferon-y by 
lymphocytes and dendritic cells [12,13]. Keratinocytes play an important role in the regulation of skin 
inflammation resulting from environmental and immune cells stimuli. Mitotic activity of basal 
keratinocytes is increased by a factor of almost 50 in psoriatic skin compared to normal skin, so that 
keratinocytes only need 3 to 5 days to move from the basal layer to the comifled layer instead of the 
normal 28 to 30 days [14]. Histological abnormalities in psoriasis have been well documented and 
include multiple elements that indicate immune-mediated inflammation. However, the cellular and 
molecular mechanisms underlying the pathogenesis of this alteration in epidermal growth and 
differentiation remains incompletely understood. 

The fundamental roles of most of the Wnt proteins in controlling cell proliferation and differentiation 
are well established. However, limited information is available about the role of Wnt signaling in 
hyper-proliferative skin diseases. Nonetheless, global gene expression analysis of psoriasis revealed 
significant changes in several members of the Wnt ligand family. Notably, non-canonical Wnts like 
WntSa, Wntl6B, and Wnt 1 OA were shown to be up-regulated. Conversely, gene expression profiling 
of canonical Wnt revealed down-regulation of nearly all members of the canonical Wnt signaling 
pathway during psoriasis [15]. Of the non-canonical Wnt proteins, Wnt5a was described to have a 
suppressive effect on normal human keratinocytes growth, while Wnt 16b was shown to induce cell 
proliferation and prolong clonogenicity in primary keratinocytes [15,16]. However, the pathogenic role 
of canonical Wnt proteins in proliferation and differentiation of keratinocytes is not completely known. 

In order to better understand the involvement of canonical Wnt signaling in proliferation and 
differentiation of keratinocytes, in this study, we performed both experimental research and a series of 
computational proteomic analyses. The effect of treatment of human primary keratinocytes with the 
canonical Wnt ligand, Wnt3a, on cell proliferation and differentiation was analyzed. Finally, 
we used computational methods to determine the relationship between the Wnt family proteins and 
their conserved regions, the aliphatic index and grand average of hydrophobicity (GRAVY) values, 
evolution origin, and protein-protein interactions. We corroborated our computational analysis results 
experimentally by looking at the proliferation and differentiation of keratinocytes. Finally, we concluded 
that Wnt3a has a significant inhibitory role in the proliferation and differentiation process of 
keratinocytes and a hypothesis was purported to explain this inhibitory role of Wnt3a in the 
proliferation and differentiation process of keratinocytes. 

2. Results 

2.1. Effect ofWntSa on the Proliferation of Cultured Keratinocytes 

Initially, we attempted to determine whether potentiating canonical Wnt signaling may have any 
effect on proliferation of human keratinocytes. Treatment of recombinant Wnt3a at a concentration of 
50 and 100 ng/mL showed inhibition of keratinocytes proliferation (Figure lA). The inhibition of 
proliferating cells was greater at the higher concentration (100 ng/mL) as compared to the lower 
concentration (50 ng/mL). The anti-proliferative effect of Wnt3a was similarly observed following 
1, 2, and 3 days of exposure. Since, expression of the human Ki-67 protein is strictly associated with 
cell proliferation [17], the anti-proliferative effect of Wnt3a on normal keratinocytes was determined 
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by measuring the expression level of Ki-67 using real-time PCR. The Ki-67 mRNA expression levels 
confirmed the inhibitory effect of Wnt3a on the proliferation of keratinocytes (Figure IB). 

Figure 1. The effect of Wnt3a on the proliferation of keratinocytes. (A) Wnt3a suppressed the 
proliferation of normal keratinocytes in a dose-dependent manner; (B) The anti-proliferative 
effect of WntSa on normal keratinocytes was confirmed using real-time PCR of Ki-67. 
'^p< 0.05, ^^p < 0.01, < 0.001 v^. control. 




2.2. Effect of Exogenous Wnt3a on the Differentiation of Cultured Keratinocytes 



In order to further analyze the effect of recombinant Wnt3a on the differentiation process of 
keratinocytes, total RNA of cultured keratinocytes was isolated after 24 h of treatment with recombinant 
WntSa (50 ng/mL). The mRNA levels of early and late differentiation markers for keratinocytes, 
including involucrin, keratin 1, loricrin, and keratin 10, were analyzed via real-time PCR. However, 
no significant change in mRNA levels were observed for any of the markers (Figure 2). 

Figure 2. Effect of exogenous Wnt3a on the differentiation process of keratinocytes. 
Quantitative real-time PCR analysis of early and late differentiation markers for 
keratinocytes (A) involucrin; (B) loricrin; (C) keratin 1; and (D) keratin 10, respectively, 
showed no significant effect of WntSa on the differentiation process of keratinocytes. 
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2.3. Effect ofWntSa on Proliferation ofTNFprocess ofkeratinocytes.a 

TNFa is a potent inflammatory cytokine that is highly expressed in psoriatic skin. Moreover, 
remarkable improvements were observed in clinical trials with TNFa antagonists [18,19] suggesting a 
crucial role of TNFa in psoriatic pathogenesis. Therefore, to determine the effect of Wnt3a on 
keratinocyte proliferation under stimulation with TNFa, human primary keratinocytes were treated 
with TNFa (2 ng/mL) in the absence or presence of Wnt3a (50 ng/mL) (Figure 3). TNFa induced 
inhibition of proliferating in normal keratinocytes for the observed time period of 1, 2 and 3 days 
(Figure 3 A). Co-treatment with TNFa and Wnt3a significantly suppressed the proliferative ability of 
keratinocytes as compared to TNFa alone. Furthermore, mRNA expression levels of Ki-67 confirmed 
this inhibitory effect of TNFa and Wnt3a on keratinocytes (Figure 3B). 

Figure 3. The effect of Wnt3a on the proUferation of TNFa-stimulated keratinocytes. (A) Wnt3a 
suppressed the proliferation of TNFa-stimulated keratinocytes; (B) The anti-proliferative 
effect of Wnt3a on TNFa-stimulated keratinocytes was confirmed using real-time PCR of 
Ki-67. ''p< 0.05, < 0.01, < o.OOl v^. control. 
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2.4. Effect ofWntSa on Differentiation of TNFa-Stimulated Keratinocytes 

In order to evaluate the effect of Wnt3a on the differentiation process of TNFa-stimulated 
keratinocytes, primary keratinocytes were treated with TNFa (2 ng/mL) with or without Wnt3a for 
24 h. Total RNA was collected and mRNA levels of involucrin, keratin 1, loricrin, and keratin 10 were 
analyzed via real-time PCR. The mRNA expression levels of early and late differentiation markers, 
keratin 1 and loricrin, showed a 2.5-3.0 fold increase after co-treatment with TNFa and Wnt3a 
(Figure 4B,C). However, no significant effect of Wnt3a was observed on the expression levels of 
involucrin and Keratin 10 (Figure 4A,D). 

2.5. Data Collection of Wnt Family Proteins for Computational Analysis 

Table SI shows the protein information related to the Wnt family proteins analysed in this study. 
The protein ID, accession, and version were obtained from the NCBI database. The 19 member 
proteins of the Wnt family consist of Wntl, Wnt2, Wnt2b, Wnt3, Wnt3a, Wnt4, Wnt5a, Wnt5b, Wnt6, 
Wnt7a, Wnt7b, WntSa, WntSb, Wnt9a, Wnt9b, WntlOa, WntlOb, Wntl 1, and Wntl6. Figure 5 depicts the 
number of amino acids present in each protein of the protein family. The highest number of amino acids 
were determined in WntlOa (417 aa), whilst the lowest number was found in Wnt6 (338 aa). Specifically, 
the number of amino acids in Wnt3a and Wnt5a were 352 and 365 amino acids, respectively. 
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Figure 4. Effect of Wnt3a on the differentiation of TNFa-stimulated keratinocytes. Early 
differentiation markers, Involucrin (A) and keratin 10 (D) showed no significant changes 
under the combinatory treatment of TNFa and Wnt3a, while mRNA levels of the early 
differentiation marker Keratin 1 (C) and late differentiation marker loricrin (B) were 
significantly increased by Wnt3a in TNF stimulated keratinocytes. p< 0.01 v^. control. 
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Figure 5. Number of amino acids in the sequence of the various Wnt family member 
proteins and the trend in their amino acid number. The proteins, WntSa and WntSa, 
are marked with an arrow. 



CI 
O 

s 

CI 



O 



4S0 
400 

tso 

2S0 
200 
ISO 
100 
SO 



» i i ♦ 



Y-0.6)(+355 
ft* = 0,0373 



I I I I I I I 



I I I I I I 



Wnt family proteins 

2.6. Multiple Sequences Alignment (MSA) of Wnt Proteins 



Clustal Omega output was visualised through JalView and the results are shown in Figure 6. 
In total, 16 aligned blocks, including 41, 98-118, 124-130, 132-138, 152-200, 245-270, 282-298, 
306-322, 322-346, 380-398, 405-423, 428-433, 436-442, 445-446, and 471-476 were found. 
The first block is the smallest containing a single residue at the 41st position, whereas the largest 
aligned block contained 48 residues (152-200). We found that Arg (R), Try (W), Asn (N), Cys (C), 
Glu (E), Gly (G), Lys (K), His (H), and Ser (S) were well aligned at different positions. 
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Figure 6. Multiple sequence alignment of the Wnt family member proteins. Wnt3a and 
Wnt5a align sequence have been marked (arrows) to show the similarity of the sequence. 
Color codes used for amino acid residues are as follows: A,I,L,M,F,W,V = Blue; 
R,K = M; N,Q ,S,T = Green; C = Blue (Threshhold + 60%); C 
E,D = Magen^; G = Orange; H,Y = ; P = Yellow. 
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2. 7. Analysis of the Aliphatic Index and GRA VY of the Wnt Family Proteins 

The aliphatic side chains are essential for the biological structure and function of proteins, and 
contain residues, such as alanine, valine, leucine, isoleucine etc. The side chains are used to describe 
the relative volume, called the aliphatic index. Aliphatic hydrophobicity is augmented with an increase 
in temperature and it is consequently a positive factor to enhance the thermal stability of globular 
proteins [20]. Our analysis showed the highest aliphatic index values in Wnt 11 (73.59). Conversely, 
the lowest value was noted in Wnt2 (63.33) (Figure 7). In this aliphatic index value scale, a low aliphatic 
index value was noted for Wnt3a (62.93) and a moderate aliphatic index value for WntSa (65.48). 

Figure 7. Aliphatic index values of Wnt family member proteins. Wnt3a and WntSa 
aliphatic index have been marked (arrows) to show the resemblance between the values. 

75 -1 
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The grand average of hydrophobicity (GRAVY) is linked with protein solubility. It has been 
observed that positive GRAVY values are positively associated with hydrophobicity and negatively 
associated with the hydrophilicity. Hydrophilic proteins form a higher number of hydrogen bonds with 
water molecules, thereby enhancing the solubility of the protein. The ProtParam GRAVY study 
predicted that all the Wnt family proteins have negative values (Figure 8A). Our study revealed that 
the highest negative GRAVY value could be established for Wnt 10a (-0.455). Conversely, the lowest 
negative GRAVY value was determined for Wnt9b (-0.265). After the segment determination of 
the GRAVY value, we found both hydrophobic and hydrophilic segments in Wnt3a and Wnt5a 
(Figure 8 A). This analysis showed that the 1st segment (1 to 70 aa) is hydrophobic in Wnt3a and 
Wnt5a. Conversely, the 3rd (141-210), 4th (121-280), and last (281-252 or 281-265) segments are 
hydrophilic in both Wnt3a and Wnt5a. 

2.8. Molecular Evolution Analysis of the Wnt Family Proteins 

To understand the molecular evolution of the Wnt family proteins, a phylogenetic tree of the Wnt 
family proteins was generated and is depicted in Figure 9. The tree illustrates bootstrap values at the 
inner nodes. It has been established that 5 different groups of proteins have 100% bootstrap 
replications. We found 100% bootstrap replications between the proteins Wnt2 and Wnt2b; Wnt9a and 
Wnt9b; Wnt5a and Wnt5b; Wnt8a and Wnt8b; Wnt3 and Wnt3a. The results also showed a 99% 
bootstrap replication value between Wnt7a and Wnt7b, as well as Wnt 10a and Wnt 10b. 
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2.9. Creation of the Phylogenetic Tree for the Wnt Family Seed Alignment 

We also constructed the "phylogenetic tree for family seed alignment" for the Wnt family proteins. 
It uses sequence similarity across species of Wnt family proteins and portrays the phylogenetic 
tree for family seed alignment, as shown in Figure 10. The tree shows sequence similarity within 
115 different proteins. 

Figure 8. Grand average of hydrophobicity (GRAVY) values of the Wnt family proteins. 
(A) The GRAVY values of 19 proteins of the Wnt family show that all the members of the 
family are hydrophilic in nature; (B) Segmental analysis of GRAVY values of two proteins: 
Wnt3a and WntSa. The analysis shows that the 1st segment (1 to 70 aa) is hydrophobic and the 
3rd (141-210), 4th (121-280), and last (281-252 or 281-265) segments are hydrophilic for 
both proteins. Arrow shows the similarity of the segments between the two proteins. 
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Figure 9. Phylogenetic tree of the 19 members of the Wnt family. WntSa and WntSa are 
marked with arrows to show their position. 
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2.10. Conservation Assessment of Wnt Family Proteins in Different Species 

The conservation of Wnt family proteins in different species was determined and is documented in 
Figure 11. These proteins are more conserved in members of the phylum Chordata. Other than the 
phylum Chordata, these proteins are also more or less conserved in other families of eukaryotes. 
The family proteins are not conserved in prokaryotes. 

Figure 10. The phylogenetic tree for Wnt family seed alignment utilizing a neighbor join 
tree algorithm with a local bootstrap based on 100 resamples (which are shown next to 
the tree nodes). 
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Figure 11. Conservation of Wnt family proteins in different species. The species range 
includes both eukaryotes and prokaryotes. Conservation of the sequences has been marked 
with different red color (from lighter red to darker red). For example, lighter red means no 
detectable homology while darker red means 100% sequence conservation. 
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2.77. Prediction of Precise Version of Sequence Similarity and Generation of HMM Logos of Wnt 
Family Proteins 

Two types of logos were created using two different tools, i.e., the sequence logo using WebLogo 
software and HMM logos using Pfam database. The generated sequence logos are shown in 
Figure 12A,B. We used 70 amino acids each time to determine the conserved part of the Wnt family 
proteins. From the sequence logo, we noted that the first sequence logo is the tallest sequence logo 
(3.4 bits). Methionine (M) may be a highly conserved amino acid in Wnt family proteins at the first 
position. In addition, we found conserved blocks at 12-17; 18-29; 32-37; 46-55; 60-64; 85-91; 
93-96; 99-107; 111-134; 135-139 (Figure 12A). From the HMM logos, we found high sequence 
logos at various positions, including the 37th (C); 46th (W); 48th (C); 68th (E); 71st (F,Y); 
88th (C) (Figure 12B). 
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Figure 12. Generated sequence logo of Wnt family proteins. (A) Sequence logo up to 
140 sequences; (B) HMM logo of Wnt family proteins up to 89 sequences. Default colour 
of webLogo and HMM logo servers has been used. 
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2. 12. Understanding the Protein-Protein Interaction Network of Wnt Family Proteins 

The protein-protein network between various Wnt family proteins is depicted in Figure 13. 
The developed three types of the protein-protein network are shown in Figure 13A (evidence 
network), Figure 13B (actions network), and Figure 13C (confidence network). The input file for the 
Wnt family proteins is shown in Figure SI. This figure shows a complex and highly interconnected 
network between the various Wnt family proteins. From the evidence network, we found that all 19 
proteins are closely connected. From the actions network, we determined that binding between 
proteins is located in the upper position, such as for Wnt3a, WntSa, Wntl, Wnt3, and Wnt2. The 
network landscape explicates that these proteins are not only interlinked among themselves, but are 
also interlinked with many other proteins, such as FZDl, FZD4, FZD5, FZDIO, WIFl, DVLl, DVL3, 
LRP5, and LRP6. The confidence network shows that all the proteins are strongly related. An overall 
stronger association was found between Wnt3a and WntSa. 

Figure 13. Protein-Protein interaction network of Wnt family proteins using STRING. 
Arrows shows the similarity in the network between the nodes. (A) Network type-evidence 
network; (B) Network type-actions network; (C) Network type-confidence network. 
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3. Discussion 

Various skin diseases, especially psoriasis, are characterized by alterations in keratinocyte growth 
and differentiation [21-23]. In normal skin, the fraction of proliferating keratinocytes is probably 
around 20%, whereas during psoriasis it is almost 100%. In addition, the mean cell cycle time is reduced 
from 13 days to 36 h [14]. These changes are associated with altered expression of genes in the 
epidermis. The late keratinocyte differentiation markers, such as filaggrin and loricrin, are down-regulated, 
while early differentiation markers, such as involucrin, are up-regulated. Keratinocytes activated by 
various stressors, including cytokines and growth factors, display a different phenotype (expressing 
K6, K16 and K17) compared to differentiated cells (expressing Kl and KIO) [11]. Recently, TNFa was 
found to modulate the expression of these proteins via the c-Jun A^-terminal kinases (JNKs) dependent 
pathway. It was reported that clinical treatment of psoriasis patients with TNFa antagonists resulted in 
significant improvement of the epidermal barrier protein expression [14]. As mentioned above, the 
Wnt family of signaling proteins is a set of highly conserved molecules that participate and control 
processes, such as cell proliferation, cell-fate differentiation, and differentiation during adult 
homeostasis [3]. To date, 19 Wnt proteins and 10 Frizzled transmembrane receptors have been 
documented. Although some knowledge has been accumulated about the state and/or function of Wnt 
signaling in hyperproliferative skin disorders, very few studies were performed to determine whether 
and how the Wnt pathway affects keratinocyte proliferation and differentiation [15,24-27]. 
A recent study, based on microarray interrogation of a large probe set and a large cohort, confirmed 
that WntSa is significantly upregulated in psoriatic lesions [15]. Expression levels of Wnt 1 OA and 
Wnt7A were found to be moderately increased, while other Wnt members were either not changed or 
showed decreased expression. Moreover, the study reported a growth-suppressive effect of WntSa on 
human keratinocytes. The effect of the Wnt canonical pathway in keratinocyte biology has been 
controversial. A recent study showed that R-spondin 2 (Rspo2), a Wnt/p-catenin signaling agonist, 
is able to increase normal keratinocyte proliferation in a dose-dependent fashion. Rspo2 was found to 
synergize with Wnt3a to bring about a significant improvement in proliferation [24]. However, another 
study was unable to observe any effect of recombinant Wnt3a on normal keratinocyte growth and 
migration [15]. Thus, based on these contradictory reports, it remains unclear whether the canonical 
pathway has any effect on the proliferation and differentiation of normal keratinocytes. Therefore, 
in this study, we investigated the effect of exogenous recombinant Wnt3a on the proliferation and 
differentiation ability of human keratinocytes. Our results indicate that treatment with Wnt3a 
suppresses the proliferative ability of human keratinocytes (Figure 1 A,B), while no effect of Wnt3a 
treatment was observed on the expression levels of early and late differentiation markers (Figure 2A-D). 

TNFa is a potent inflammatory cytokine that is highly expressed in psoriatic skin. This cytokine has 
a crucial role in psoriasis pathogenesis, as demonstrated by the efficacy of TNFa-targeted therapeutics. 
TNFa is a powerful inducer of inflammatory gene products in human keratinocytes [28,29]. 
To investigate the effect of Wnt3a on keratinocytes in psoriatic skin lesions, we induced inflammatory 
conditions in keratinocytes by treating them with TNFa. TNF receptor signaling has been shown to 
induce psoriasis-like skin inflammation in keratinocytes [30]. Under these conditions, Wnt3a was able 
to suppress proliferation of keratinocytes (Figure 3A,B). Taken together, our study results show a 
suppressive role of Wnt3a on the growth of normal and TNFa-stimulated keratinocytes. Therefore, 
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under the conditions that exist in psoriatic lesions, application of Wnt3a might be applied as an 
anti-proliferative agent to control the growth of abnormal keratinocytes. 

Moreover, Wnt3a treatment increased the expression levels of early and late differentiation markers 
in TNFa-stimulated keratinocytes, such as keratin 1 and loricrin (Figure 4B,C). Though the expression 
levels of involucrin and keratin 10 were found to be unaltered, increased expression levels of keratin 1 
and loricrin induced by Wnt3a suggest initiation of the differentiation process in TNFa-stimulated 
keratinocytes. These results imply that WntSa may have a recovering effect on the altered differentiation 
process in keratinocytes, as observed in psoriatic skin lesions. 

Our in silico analysis showed that the whole Wnt family of signaling proteins are highly conserved. 
Due to this conserved nature, Wnt proteins are closely associated with vital cellular processes, such as 
cell signaling, cell-cell communication, and cell fate determination [31,32]. These highly conserved 
proteins are not only found in humans, but also in other species. Our results indicate that the Wnt 
family is highly conserved among several members of the phylum Chordata, which resembles the Wnt 
signal-transduction pathway which is widely conserved and may be necessary for important biological 
process (Figures 10 and 11). The Wnt signaling pathway has not only been found to be conserved in 
Chordata, but even in other species, such as C. elegans, Drosophila etc. [33,34]. Furthermore, amino acid 
comparison of the 19 protein members of the Wnt family has shown a homologues range of similarity 
from 27% to 83% [35]. In our analysis, we also observed high sequence similarity within all Wnt 
proteins (Figure 6). Sequence similarity between Wnt3a and WntSa was also observed in several 
conserved regions. Conserved residues may be associated with the formation of U-shaped binding 
contours, which may be necessarily needed to act as a signaling cascade and help in receptor binding [36]. 
The growth-suppressive effect of WntSa on normal human keratinocytes has already been reported [15]. 
As demonstrated by our in silico analysis, a similarity between structural and functional properties 
of Wnt3a and WntSa was observed. Comparable resemblances were observed in the amino acids 
sequence (Figure 5), multiple sequence alignment (Figure 6), aliphatic index (Figure 7), and GRAVY, 
as well as segmental GRAVY values (Figure 8). Moreover, we found a strong association between 
Wnt3a and WntSa from the protein-protein interaction networks (evidence network, actions network, 
and confidence network) during different biological processes (Figure 13). This strongly demonstrates 
the existence of mutually connected networks between these two proteins. Taken together, the results 
obtained from in silico analysis shows a close relationship between Wnt3a and WntSa, and therefore 
supports our experimental result that Wnt3a is indeed associated with proliferation and differentiation 
of keratinocytes along with WntSa. However, more experimental studies are necessary to fully 
understand the association between Wnt3a and WntSa, as well as to determine the role of the Wnt3a 
signaling cascade in keratinocyte proliferation and differentiation. 

Wnt proteins are not only highly conserved, but their aliphatic index also shows that they are highly 
stable (Figure 7). Our GRAVY value shows that these proteins are hydrophilic in nature. However, 
several hydrophobic patches were also detected in these proteins (Figure 8A). The first two segments 
in Wnt3a were found to be hydrophobic, while the other three were hydrophilic in nature (Figure 8B). 
Finch et al. [37] stated that the signal peptide part of the NH2-terminus of a secreted protein contains 
hydrophobic amino acid segments, which are usually 27 amino acids in number. This observation is 
quite similar to the pattern of hydrophobic amino acids observed by us, indicating that the first 
segment of Wnt3a might be acting as a signaling peptide for the protein. 



Int. J. Mol Set 2014, 15 



5487 



One US patent on human keratinocytes supports the fact that keratinocytes are maintained by 
hydrophilic membrane [38]. Similar observations on the hydrophihc character of the membrane of 
human keratinocytes were also reported by Yamato et al, [39]. On the basis of the above findings, it may 
be speculated that the hydrophobic segment present in Wnt3a may help its attachment to the 
hydrophilic membrane of keratinocytes. It is well established that non-specific hydrophobic 
associations exist within the plasma membrane. For example, hydrophobic and non-hydrophobic 
interactions (electrostatic attractions) were reported as a significant event in the interaction of the 
histactophilin protein with natural hydrophobic anchors [40]. Alternatively, a possible of interaction of 
the hydrophobic segment of Wnt, acting as signal peptide, with hydrophobic portion (lipid bilayer) of 
the membrane may be purported. Studies by Gonnet et al. have identified hydrophobic proteins 
attached to cytoplasmic membranes of eukaryotic keratinocytes in complex biological samples [41]. 
Therefore, we developed a hypothetical structure for illustrating the interaction between the 
hydrophobic parts of Wnt3a and the hydrophilic membrane of keratinocytes (Figure 14). Here, we tried to 
depict a proposed mechanism of action to describe how Wnt3a might inhibit the process of proliferation 
and may contribute to differentiation of keratinocytes. Due to hydrophobic and hydrophilic 
interactions, Wnt3a proteins might be able to interact with the membrane of keratinocytes making 
aggregates on its surface and thus creating a hydrophilic environment outside the cell. Previously, 
hydrophilic repulsion between two biopolymers has been reported through both in vitro and in vivo 
experiments [42,43]. Therefore, it is expected that two keratinocytes with aggregated Wnt3a protein on 
their surfaces would repel each other and thus may contribute to the biological function of keratinocytes, 
as observed in the experimental results (Figure 13). Nevertheless, our hypothesis is a probable one, but is 
only possible if the hydrophilic domain contains a net positive or negative charge; the energy needed 
for repulsion. So, further studies are required to ascertain and understand the biophysical and 
biochemical interactions between these signaling molecules. Apart from this hypothetical model, 
activation of the Wnt signaling cascade may also account for suppressed proliferation, both in TNFa 
treated or untreated keratinocytes, and promote differentiation in TNFa stimulated-keratinocytes. 
Consequently, to confirm our results and proposed hypothesis and to unfold the pathway of Wnt3a, 
more studies are needed. These studies should aim to answer different questions, such as (1) whether 
Wnt3a shows similar effects in psoriatic lesions; (2) determine detailed mechanisms underlying the 
anti-proliferative effect of Wnt3a on keratinocytes; and (3) elucidate the exact role of the Wnt pathway 
in the pathogenesis of psoriasis. 

4. Methods 

4.1. Tissue Specimens 

Primary normal human skin samples were obtained from consented patients who underwent 
elective surgery at Hallym University Chuncheon Sacred Heart Hospital (Chuncheon, Korea). Sample 
collections were obtained after patients agreed to donate their skin for research purposes (read and 
signed the informed consent forms). This study was approved by the Research Ethical Committee of 
Hallym University, College of Medicine, Chuncheon, Korea. 
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Figure 14. Proposed hypothesis for the interaction of the Wnt3a protein and hydrophihc 
membrane keratinocytes and how Wnt3a produces the suppressive effect in keratinocytes. 
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4.2. Cell Culture 



Normal human skin samples were washed in saline, minced, and then treated with dispase 
(Boehringer-Mannheim, Mannheim, Germany) for 4 h at 4 "^C. The epidermis was separated and 
placed in a solution containing 0.025% trypsin (GIBCO-BRL, Grand Island, NY, USA) and 
0.01% EDTA at 37 °C for 15 min. After vigorous pipetting, the cells were pelleted and resuspended in 
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serum-free keratinocyte growth medium (KBG gold medium, GIBCO-BRL, Grand-Island, NY, USA) 
supplemented with bovine pituitary extract, recombinant human epidermal growth factor, insulin, 
and hydrocortisone. The cells were added to 100 mm dishes (Costar, Pittsburg, NY, USA) and 
incubated at 37 °C in 5% CO2. At 70%-80% confluence, the cells were passaged. 

4.3. Cell Viability Assay (MTT Assay) 

Keratinoc3^es were cultured in 96-well plates at a concentration of 1 x lo"^ cells/well for 12 h. 
Cells were then treated with PBS containing TNF-a (2 ng/mL) or Wnt3a (50 ng/mL; Tocris, Bristol, 
UK) alone, or in combination in KBG gold medium for 24, 48 and 72 h at 37 °C in 5% CO2. 
Ten micro liters of 3-(4,5-dimethlthiazol-2-yl)-2,5-diphenyl-tetrazoliumbromide (MTT) (5 mg/mL 
dissolved in PBS; Sigma, St. Lois, MO, USA) was added to each well, and the plates were incubated at 
37 °C for 2 h. The supernatant was discarded and 200 |iL of dimethyl sulfoxide (DMSO; Sigma, 
St. Lois, MO, USA) was added to dissolve the blue insoluble MTT formazan produced by 
mitochondrial succinate dehydrogenase. The absorbance was measured spectrophotometrically at 
570 nm. Experiments were repeated three times and the data were expressed as the means =t SD. 

4.4. Real Time PCR 

The cultured keratinocytes were treated with PBS containing TNF-a (2 ng/mL) or Wnt3a (50 ng/mL) 
alone or in combination in KBG gold medium for 24 h at 37 °C in 5% CO2. Total RNA was extracted 
by using the Trizol reagent (Invitrogen, San Diego, CA, USA) and cDNA was synthesized with 2 |ig 
of total RNA using Superscript II (Invitrogen, San Diego, CA, USA) according to the manufacturer's 
instructions. Real-time PCR was performed using 1 |iL of cDNA in a 20 |iL reaction volume with the 
Rotor gene system and QuantiTect SYBR Green PCR Master Mix (Bioneer, Daejeon, Korea). 
The temperature profile of the reaction was 95 °C for 15 min, 45 cycles of denaturation at 94 °C for 
20 s, annealing at 60 "^C for 20 s, and extension at 72 °C for 25 s. Primers for keratin 1, keratin 10, 
involucrin, loricrin, and Ki-67 were obtained from Bioneer (Daejeon, Korea). The forward and reverse 
primer sequences are listed in Table 1. The relative mRNA levels were normalized by using GAPDH as a 
housekeeping gene. Rotor gene software was used to compare each gene sample level. 



Table 1. Primers for real-time PCR. 



Target 


Forward primer (5-3*) 


Reverse primer (3-5*) 


GenBank No. 


Loricrin 


TAGAGCTCTCATGATGCTACCCGA 


GGGTTGGGAGGTAGTTGTACAGAA 


NM_000427.2 


Involucrin 


AATGAAACAGCCAACTCCACTGCC 


AGGCAGTCCCTTTACAGCAGTCAT 


NM_005517.2 


Keratin 1 


GGCAGTTCCAGCGTGAAGTTTGTT 


TTCTCCGGTAAGGCTGGGACAAAT 


NM_006121.3 


Keratin 10 


TTGGTGGAGGTAGCTTTCGTGGAA 


AGAAGGCCACCATCTCCTCCAAAT 


NM_000421.3 


Ki-67 


TGACAAGCCCACGACTGATGAGAA 


CTTTGCCTGCTGATGGTGTTCGTT 


NM_002417.4 


GAPDH 


TCGACAGTCAGCCGCATCTTCTTT 


ACCAAATCCGTTGACTCCGACCTT 


NM_002046.3 



4.5. Data Collection of Wnt Family Proteins for Computational Analysis 

Computational analysis of Wnt family proteins was performed by obtaining sequence information 
of human Wnt family proteins from the NCBI database (http://www.ncbi.nlm.nih.gov) [44]. The 
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functional protein sequences were collected in FASTA format along with the accession number from 
the NCBI database for further analysis. 

4.6. Multiple Sequences Alignment (MSA) of Wnt Proteins 

To understand the resemblance within Wnt proteins, Clustal Omega was used. Clustal Omega 
server has a graphical edge [45] and utilises several other algorithms, such as "progressive algorithm" 
and Hidden Markov Models (HMMs) profile to generate alignments. In fact, HMM, is the formal basis 
for the creation of probabilistic models of linear sequence labelling problems [46,47] and also uses pair 
wise comparison methods for large-scale sequence analysis. The graphical yield of MSA was 
visualised through JalView (http://www.jalview.org/). 

4. 7. Analysis of Aliphatic Index and GRA VY of Wnt Family Proteins 

A comparison of the various physical and chemical parameters of Wnt family proteins was 
performed. Two physicochemical parameters, i.e., the aliphatic index and grand average of 
hydrophobicity (GRAVY), were calculated using the ProtParam software from the ExPASy server [48]. 
In order to understand the segment related hydrophobicity, hydrophilicity and similarity of the two 
proteins (Wnt3a and WntSa), GRAVY values of different segments of the protein chain were 
determined. Each segment selected for analysis contained 70 amino acids (aa) (1st segment, 1 to 70 aa), 
(2nd segment, 70 to 140 aa), (3rd segment, 141 to 210 aa), (4th segment, 121 to 280 aa) and last 
(5th segment, 281-252 or 281-265 aa) of the WntSa and Wnt5a proteins. 

4.8. Molecular Evolution Analysis of Wnt Family Proteins 

The molecular evolution was analyzed by developing a phylogenetic tree using the "Phylogeny.fr" 
server. This server uses several significant algorithms, such as MUSCLE (https://www.ebi.ac.uk/ 
Tools/msa/muscle/) for multiple alignments, PhyML (http://atgc.lirmm.fr/phyml) for tree building, and 
TreeDyn (http://www.treedyn.org) for tree rendering, as well as G blocks to reconstruct a robust 
phylogenetic tree from a set of sequences [49]. This also uses the maximum-likelihood method for tree 
construction along with a Bayesian method for molecular phylogenetics. 

4.9. Creation of the Phylogenetic Tree for Wnt Family Seed Alignment 

To understand the distribution of Wnt family proteins across different species, Pfam database was 
used [50]. This database contains a large collection of protein families and domain databases, which 
facilitates the development of a phylogenetic tree for family seed alignment. It utilizes FastTree 
software, which in turn uses approximately-maximum-likelihood phylogenetic trees. Neighbor join 
trees were developed using local bootstrap based on 100 resamples to calculate the phylogenetic tree 
for family seed alignment. 

4.10. Conservation Assessment of Wnt Family Proteins in Different Species 

STRING server was used to understand the conservation of Wnt family proteins in different 
species [51]. STRING systematically carries out orthology transfers, utilizing two systems at a time: 
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precomputed orthologs from the COG database [52] and a homology-based orthology scheme, 
computed de novo [53]. This server uses more than 7x10^^ pair wise protein comparisons, using the 
sensitive Smith- Waterman dynamic programming algorithm [54]. 

4.11. Prediction of Precise Version of Sequence Similarity through Sequence Logo and Generation of 
HMM Logos of Wnt Family Proteins 

To understand the more affluent and accurate description of sequence similarity, WebLogo 
server [55] was used. This tool allows the graphic representation of amino acids and also highlights the 
patterns in a set of aligned sequences [55]. It has been extensively documented that the sequence logo 
is at an exact position in the alignment of the residues [56], which can be described as i?seq according to: 

N 

^seq = ^max-^obs= ^^^2^ ~ C^P^^Og^p^) (1) 

where the difference 7?seq is between the maximum possible entropy and the entropy of the observed 
symbol distribution, /?n is the observed frequency of symbol w at a particular sequence position, and 
N diXQ the number of distinct symbols for the given sequence type. 

Presently, the Profile Hidden Markov Model (pHMMs) is a widely used tool for protein family 
research. In this work, for the generation of HMM logos, we used Pfam database [50], which provides 
a quick overview of the properties of a HMM in a graphical format [57]. 

4. 12. Understanding the Protein-Protein Interaction Network of Wnt Family Proteins 

Three types of protein-protein interaction network were developed using the STRING server for 
the 19 Wnt family proteins: confidence network, evidence network, and actions network. This web-based 
server investigates possible protein-protein interactions by utilizing direct (physical) and indirect 
(functional) associations among the members [51]. 

4.13. Statistical Analysis 

All the statistical data were analyzed via Graphpad Prism 5.0 (GraphPad Software, San Diego, CA, 
USA) and evaluated by the two-tailed Student t test. Values of < 0.05 were considered to indicate 
statistical significance. 

5. Conclusions 

Here in, we observed that incubation of human keratinocytes with the canonical Wnt protein, 
Wnt3a, suppresses the proliferative ability of keratinocytes, while no such conclusive effect was 
observed on the expression level of differentiation markers. Addition of Wnt3a to TNFa-stimulated 
keratinocytes also resulted in an anti-proliferative response in keratinocytes. Wnt3a incubation under 
TNFa-stimulated conditions was able to induce the expression of early and late differentiation 
markers, implicating WntSa in the regulation of the differentiation process under pathogenic 
conditions. Computational analysis of the Wnt ligand family highlighted biochemical similarities 
between WntSa and Wnt5a. The non-canonical Wnt member, Wnt5a, is a well-known suppressor of 



Int. J. Mol Set 2014, 15 



5492 



keratinocytes proliferation. Moreover, on the basis of the GRAVY value, we were able to deduce a 
hypothetical model for WntSa's anti-proliferative effect on keratinocytes. Our results and proposed 
hypothesis might be useful for future researchers to obtain fundamental information about the 
structural and functional significance relating to the evolutionary biology of Wnt proteins and their 
role in various disease states, such as psoriasis. 
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